\(\int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx\) [58]

Optimal result

Integrand size = 28, antiderivative size = 177 \[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\frac {2 a^{5/2} c^2 \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{f}-\frac {2 a^3 c^2 \tan (e+f x)}{f \sqrt {a+a \sec (e+f x)}}+\frac {2 a^4 c^2 \tan ^3(e+f x)}{3 f (a+a \sec (e+f x))^{3/2}}+\frac {6 a^5 c^2 \tan ^5(e+f x)}{5 f (a+a \sec (e+f x))^{5/2}}+\frac {2 a^6 c^2 \tan ^7(e+f x)}{7 f (a+a \sec (e+f x))^{7/2}} \]

[Out]

2*a^(5/2)*c^2*arctan(a^(1/2)*tan(f*x+e)/(a+a*sec(f*x+e))^(1/2))/f-2*a^3*c^2*tan(f*x+e)/f/(a+a*sec(f*x+e))^(1/2
)+2/3*a^4*c^2*tan(f*x+e)^3/f/(a+a*sec(f*x+e))^(3/2)+6/5*a^5*c^2*tan(f*x+e)^5/f/(a+a*sec(f*x+e))^(5/2)+2/7*a^6*
c^2*tan(f*x+e)^7/f/(a+a*sec(f*x+e))^(7/2)

Rubi [A] (verified)

Time = 0.22 (sec) , antiderivative size = 177, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {3989, 3972, 472, 209} \[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\frac {2 a^{5/2} c^2 \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {a \sec (e+f x)+a}}\right )}{f}+\frac {2 a^6 c^2 \tan ^7(e+f x)}{7 f (a \sec (e+f x)+a)^{7/2}}+\frac {6 a^5 c^2 \tan ^5(e+f x)}{5 f (a \sec (e+f x)+a)^{5/2}}+\frac {2 a^4 c^2 \tan ^3(e+f x)}{3 f (a \sec (e+f x)+a)^{3/2}}-\frac {2 a^3 c^2 \tan (e+f x)}{f \sqrt {a \sec (e+f x)+a}} \]

[In]

Int[(a + a*Sec[e + f*x])^(5/2)*(c - c*Sec[e + f*x])^2,x]

[Out]

(2*a^(5/2)*c^2*ArcTan[(Sqrt[a]*Tan[e + f*x])/Sqrt[a + a*Sec[e + f*x]]])/f - (2*a^3*c^2*Tan[e + f*x])/(f*Sqrt[a
+ a*Sec[e + f*x]]) + (2*a^4*c^2*Tan[e + f*x]^3)/(3*f*(a + a*Sec[e + f*x])^(3/2)) + (6*a^5*c^2*Tan[e + f*x]^5)
/(5*f*(a + a*Sec[e + f*x])^(5/2)) + (2*a^6*c^2*Tan[e + f*x]^7)/(7*f*(a + a*Sec[e + f*x])^(7/2))

Rule 209

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*ArcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /;
FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])

Rule 472

Int[(((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_))/((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Int[ExpandIntegr
and[(e*x)^m*((a + b*x^n)^p/(c + d*x^n)), x], x] /; FreeQ[{a, b, c, d, e, m}, x] && NeQ[b*c - a*d, 0] && IGtQ[n
, 0] && IGtQ[p, 0] && (IntegerQ[m] || IGtQ[2*(m + 1), 0] ||  !RationalQ[m])

Rule 3972

Int[cot[(c_.) + (d_.)*(x_)]^(m_.)*(csc[(c_.) + (d_.)*(x_)]*(b_.) + (a_))^(n_.), x_Symbol] :> Dist[-2*(a^(m/2 +
n + 1/2)/d), Subst[Int[x^m*((2 + a*x^2)^(m/2 + n - 1/2)/(1 + a*x^2)), x], x, Cot[c + d*x]/Sqrt[a + b*Csc[c +
d*x]]], x] /; FreeQ[{a, b, c, d}, x] && EqQ[a^2 - b^2, 0] && IntegerQ[m/2] && IntegerQ[n - 1/2]

Rule 3989

Int[(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_.)*(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_))^(n_.), x_Symbol] :> Di
st[((-a)*c)^m, Int[Cot[e + f*x]^(2*m)*(c + d*Csc[e + f*x])^(n - m), x], x] /; FreeQ[{a, b, c, d, e, f, n}, x]
&& EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && IntegerQ[m] && RationalQ[n] &&  !(IntegerQ[n] && GtQ[m - n, 0])

Rubi steps \begin{align*} \text {integral}& = \left (a^2 c^2\right ) \int \sqrt {a+a \sec (e+f x)} \tan ^4(e+f x) \, dx \\ & = -\frac {\left (2 a^5 c^2\right ) \text {Subst}\left (\int \frac {x^4 \left (2+a x^2\right )^2}{1+a x^2} \, dx,x,-\frac {\tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{f} \\ & = -\frac {\left (2 a^5 c^2\right ) \text {Subst}\left (\int \left (-\frac {1}{a^2}+\frac {x^2}{a}+3 x^4+a x^6+\frac {1}{a^2 \left (1+a x^2\right )}\right ) \, dx,x,-\frac {\tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{f} \\ & = -\frac {2 a^3 c^2 \tan (e+f x)}{f \sqrt {a+a \sec (e+f x)}}+\frac {2 a^4 c^2 \tan ^3(e+f x)}{3 f (a+a \sec (e+f x))^{3/2}}+\frac {6 a^5 c^2 \tan ^5(e+f x)}{5 f (a+a \sec (e+f x))^{5/2}}+\frac {2 a^6 c^2 \tan ^7(e+f x)}{7 f (a+a \sec (e+f x))^{7/2}}-\frac {\left (2 a^3 c^2\right ) \text {Subst}\left (\int \frac {1}{1+a x^2} \, dx,x,-\frac {\tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{f} \\ & = \frac {2 a^{5/2} c^2 \arctan \left (\frac {\sqrt {a} \tan (e+f x)}{\sqrt {a+a \sec (e+f x)}}\right )}{f}-\frac {2 a^3 c^2 \tan (e+f x)}{f \sqrt {a+a \sec (e+f x)}}+\frac {2 a^4 c^2 \tan ^3(e+f x)}{3 f (a+a \sec (e+f x))^{3/2}}+\frac {6 a^5 c^2 \tan ^5(e+f x)}{5 f (a+a \sec (e+f x))^{5/2}}+\frac {2 a^6 c^2 \tan ^7(e+f x)}{7 f (a+a \sec (e+f x))^{7/2}} \\ \end{align*}

Mathematica [A] (verified)

Time = 1.10 (sec) , antiderivative size = 124, normalized size of antiderivative = 0.70 \[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\frac {2 a^3 c^2 \left (105 \sqrt {c} \text {arctanh}\left (\frac {\sqrt {c-c \sec (e+f x)}}{\sqrt {c}}\right )+\sqrt {c-c \sec (e+f x)} \left (-92-46 \sec (e+f x)+18 \sec ^2(e+f x)+15 \sec ^3(e+f x)\right )\right ) \tan (e+f x)}{105 f \sqrt {a (1+\sec (e+f x))} \sqrt {c-c \sec (e+f x)}} \]

[In]

Integrate[(a + a*Sec[e + f*x])^(5/2)*(c - c*Sec[e + f*x])^2,x]

[Out]

(2*a^3*c^2*(105*Sqrt[c]*ArcTanh[Sqrt[c - c*Sec[e + f*x]]/Sqrt[c]] + Sqrt[c - c*Sec[e + f*x]]*(-92 - 46*Sec[e +
f*x] + 18*Sec[e + f*x]^2 + 15*Sec[e + f*x]^3))*Tan[e + f*x])/(105*f*Sqrt[a*(1 + Sec[e + f*x])]*Sqrt[c - c*Sec
[e + f*x]])

Maple [A] (verified)

Time = 22.84 (sec) , antiderivative size = 207, normalized size of antiderivative = 1.17

method	result	size
default	\(-\frac {2 a^{2} c^{2} \sqrt {a \left (\sec \left (f x +e \right )+1\right )}\, \left (-105 \,\operatorname {arctanh}\left (\frac {\sin \left (f x +e \right )}{\left (\cos \left (f x +e \right )+1\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}}\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}\, \cos \left (f x +e \right )-105 \,\operatorname {arctanh}\left (\frac {\sin \left (f x +e \right )}{\left (\cos \left (f x +e \right )+1\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}}\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}+92 \sin \left (f x +e \right )+46 \tan \left (f x +e \right )-18 \sec \left (f x +e \right ) \tan \left (f x +e \right )-15 \sec \left (f x +e \right )^{2} \tan \left (f x +e \right )\right )}{105 f \left (\cos \left (f x +e \right )+1\right )}\)	\(207\)
parts	\(\frac {2 c^{2} a^{2} \sqrt {a \left (\sec \left (f x +e \right )+1\right )}\, \left (3 \,\operatorname {arctanh}\left (\frac {\sin \left (f x +e \right )}{\left (\cos \left (f x +e \right )+1\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}}\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}\, \cos \left (f x +e \right )+3 \,\operatorname {arctanh}\left (\frac {\sin \left (f x +e \right )}{\left (\cos \left (f x +e \right )+1\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}}\right ) \sqrt {-\frac {\cos \left (f x +e \right )}{\cos \left (f x +e \right )+1}}+8 \sin \left (f x +e \right )+\tan \left (f x +e \right )\right )}{3 f \left (\cos \left (f x +e \right )+1\right )}+\frac {2 c^{2} a^{2} \left (46 \cos \left (f x +e \right )^{3}+23 \cos \left (f x +e \right )^{2}+12 \cos \left (f x +e \right )+3\right ) \sqrt {a \left (\sec \left (f x +e \right )+1\right )}\, \tan \left (f x +e \right ) \sec \left (f x +e \right )^{2}}{21 f \left (\cos \left (f x +e \right )+1\right )}-\frac {4 c^{2} a^{2} \sqrt {a \left (\sec \left (f x +e \right )+1\right )}\, \left (43 \sin \left (f x +e \right )+14 \tan \left (f x +e \right )+3 \sec \left (f x +e \right ) \tan \left (f x +e \right )\right )}{15 f \left (\cos \left (f x +e \right )+1\right )}\)	\(317\)

[In]

int((a+a*sec(f*x+e))^(5/2)*(c-c*sec(f*x+e))^2,x,method=_RETURNVERBOSE)

[Out]

-2/105*a^2*c^2/f*(a*(sec(f*x+e)+1))^(1/2)/(cos(f*x+e)+1)*(-105*arctanh(sin(f*x+e)/(cos(f*x+e)+1)/(-cos(f*x+e)/
(cos(f*x+e)+1))^(1/2))*(-cos(f*x+e)/(cos(f*x+e)+1))^(1/2)*cos(f*x+e)-105*arctanh(sin(f*x+e)/(cos(f*x+e)+1)/(-c
os(f*x+e)/(cos(f*x+e)+1))^(1/2))*(-cos(f*x+e)/(cos(f*x+e)+1))^(1/2)+92*sin(f*x+e)+46*tan(f*x+e)-18*sec(f*x+e)*
tan(f*x+e)-15*sec(f*x+e)^2*tan(f*x+e))

Fricas [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 409, normalized size of antiderivative = 2.31 \[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\left [\frac {105 \, {\left (a^{2} c^{2} \cos \left (f x + e\right )^{4} + a^{2} c^{2} \cos \left (f x + e\right )^{3}\right )} \sqrt {-a} \log \left (\frac {2 \, a \cos \left (f x + e\right )^{2} - 2 \, \sqrt {-a} \sqrt {\frac {a \cos \left (f x + e\right ) + a}{\cos \left (f x + e\right )}} \cos \left (f x + e\right ) \sin \left (f x + e\right ) + a \cos \left (f x + e\right ) - a}{\cos \left (f x + e\right ) + 1}\right ) - 2 \, {\left (92 \, a^{2} c^{2} \cos \left (f x + e\right )^{3} + 46 \, a^{2} c^{2} \cos \left (f x + e\right )^{2} - 18 \, a^{2} c^{2} \cos \left (f x + e\right ) - 15 \, a^{2} c^{2}\right )} \sqrt {\frac {a \cos \left (f x + e\right ) + a}{\cos \left (f x + e\right )}} \sin \left (f x + e\right )}{105 \, {\left (f \cos \left (f x + e\right )^{4} + f \cos \left (f x + e\right )^{3}\right )}}, -\frac {2 \, {\left (105 \, {\left (a^{2} c^{2} \cos \left (f x + e\right )^{4} + a^{2} c^{2} \cos \left (f x + e\right )^{3}\right )} \sqrt {a} \arctan \left (\frac {\sqrt {\frac {a \cos \left (f x + e\right ) + a}{\cos \left (f x + e\right )}} \cos \left (f x + e\right )}{\sqrt {a} \sin \left (f x + e\right )}\right ) + {\left (92 \, a^{2} c^{2} \cos \left (f x + e\right )^{3} + 46 \, a^{2} c^{2} \cos \left (f x + e\right )^{2} - 18 \, a^{2} c^{2} \cos \left (f x + e\right ) - 15 \, a^{2} c^{2}\right )} \sqrt {\frac {a \cos \left (f x + e\right ) + a}{\cos \left (f x + e\right )}} \sin \left (f x + e\right )\right )}}{105 \, {\left (f \cos \left (f x + e\right )^{4} + f \cos \left (f x + e\right )^{3}\right )}}\right ] \]

[In]

integrate((a+a*sec(f*x+e))^(5/2)*(c-c*sec(f*x+e))^2,x, algorithm="fricas")

[Out]

[1/105*(105*(a^2*c^2*cos(f*x + e)^4 + a^2*c^2*cos(f*x + e)^3)*sqrt(-a)*log((2*a*cos(f*x + e)^2 - 2*sqrt(-a)*sq
rt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)*sin(f*x + e) + a*cos(f*x + e) - a)/(cos(f*x + e) + 1)) - 2*
(92*a^2*c^2*cos(f*x + e)^3 + 46*a^2*c^2*cos(f*x + e)^2 - 18*a^2*c^2*cos(f*x + e) - 15*a^2*c^2)*sqrt((a*cos(f*x
+ e) + a)/cos(f*x + e))*sin(f*x + e))/(f*cos(f*x + e)^4 + f*cos(f*x + e)^3), -2/105*(105*(a^2*c^2*cos(f*x + e
)^4 + a^2*c^2*cos(f*x + e)^3)*sqrt(a)*arctan(sqrt((a*cos(f*x + e) + a)/cos(f*x + e))*cos(f*x + e)/(sqrt(a)*sin
(f*x + e))) + (92*a^2*c^2*cos(f*x + e)^3 + 46*a^2*c^2*cos(f*x + e)^2 - 18*a^2*c^2*cos(f*x + e) - 15*a^2*c^2)*s
qrt((a*cos(f*x + e) + a)/cos(f*x + e))*sin(f*x + e))/(f*cos(f*x + e)^4 + f*cos(f*x + e)^3)]

Sympy [F]

\[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=c^{2} \left (\int a^{2} \sqrt {a \sec {\left (e + f x \right )} + a}\, dx + \int \left (- 2 a^{2} \sqrt {a \sec {\left (e + f x \right )} + a} \sec ^{2}{\left (e + f x \right )}\right )\, dx + \int a^{2} \sqrt {a \sec {\left (e + f x \right )} + a} \sec ^{4}{\left (e + f x \right )}\, dx\right ) \]

[In]

integrate((a+a*sec(f*x+e))**(5/2)*(c-c*sec(f*x+e))**2,x)

[Out]

c**2*(Integral(a**2*sqrt(a*sec(e + f*x) + a), x) + Integral(-2*a**2*sqrt(a*sec(e + f*x) + a)*sec(e + f*x)**2,
x) + Integral(a**2*sqrt(a*sec(e + f*x) + a)*sec(e + f*x)**4, x))

Maxima [F]

\[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\int { {\left (a \sec \left (f x + e\right ) + a\right )}^{\frac {5}{2}} {\left (c \sec \left (f x + e\right ) - c\right )}^{2} \,d x } \]

[In]

integrate((a+a*sec(f*x+e))^(5/2)*(c-c*sec(f*x+e))^2,x, algorithm="maxima")

[Out]

-1/210*(105*((a^2*c^2*cos(2*f*x + 2*e)^2 + a^2*c^2*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*cos(2*f*x + 2*e) + a^2*c^2)*
arctan2((cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)*sin(1/2*arctan2(sin(2*f*x + 2
*e), cos(2*f*x + 2*e) + 1)), (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)*cos(1/2*
arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) + 1) - (a^2*c^2*cos(2*f*x + 2*e)^2 + a^2*c^2*sin(2*f*x + 2*e)
^2 + 2*a^2*c^2*cos(2*f*x + 2*e) + a^2*c^2)*arctan2((cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*
e) + 1)^(1/4)*sin(1/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)), (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)
^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)*cos(1/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - 1) - 2*(a^2*c^2*
f*cos(2*f*x + 2*e)^2 + a^2*c^2*f*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*f*cos(2*f*x + 2*e) + a^2*c^2*f)*integrate((((c
os(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)
*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*cos(9/2*arctan2(sin(2*f*x + 2*e)
, cos(2*f*x + 2*e))) + (cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x +
6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*sin(9/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e
))))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - ((cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*
f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*cos(9/2
*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))) - (cos(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*
f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + si
n(2*f*x + 2*e)^2)*sin(9/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(
2*f*x + 2*e) + 1)))/(((cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*
cos(2*f*x + 2*e) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1
)*cos(4*f*x + 4*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) +
1)*sin(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^
2 + (2*cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2
*e)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e)
+ 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*
f*x + 2*e)^2 + 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x +
2*e)^3 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x +
2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*
e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x
+ 2*e) + 1))^2 + (cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(
2*f*x + 2*e) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*co
s(4*f*x + 4*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*s
in(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 +
(2*cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*
sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2
*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x
+ 2*e)^2 + 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e
)^3 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e
)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2
+ 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2
*e) + 1))^2)*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)), x) + 8*(a^2*c^2*f*cos(
2*f*x + 2*e)^2 + a^2*c^2*f*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*f*cos(2*f*x + 2*e) + a^2*c^2*f)*integrate((((cos(6*f
*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2
*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*cos(7/2*arctan2(sin(2*f*x + 2*e), cos(
2*f*x + 2*e))) + (cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*s
in(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*sin(7/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*c
os(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - ((cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x +
2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*cos(7/2*arcta
n2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))) - (cos(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x +
2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*
x + 2*e)^2)*sin(7/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x
+ 2*e) + 1)))/(((cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*
f*x + 2*e) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(
4*f*x + 4*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin
(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2
*cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*si
n(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*c
os(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x +
2*e)^2 + 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^
3 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^
2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 +
2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e
) + 1))^2 + (cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x
+ 2*e) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*
x + 4*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f
*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos
(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*
f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2
*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)
^2 + 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 +
2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 +
2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*c
os(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) +
1))^2)*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)), x) - 12*(a^2*c^2*f*cos(2*f*x
+ 2*e)^2 + a^2*c^2*f*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*f*cos(2*f*x + 2*e) + a^2*c^2*f)*integrate((((cos(6*f*x +
6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x
+ 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x
+ 2*e))) + (cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*
f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*cos(5/
2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - ((cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*
sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*cos(5/2*arctan2(si
n(2*f*x + 2*e), cos(2*f*x + 2*e))) - (cos(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e)
+ cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2
*e)^2)*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e
) + 1)))/(((cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x +
2*e) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x
+ 4*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*
x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(
2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f
*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*
f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^
2 + 2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2
*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2
*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*co
s(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1
))^2 + (cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e
) + 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4
*e)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*x +
6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(2*f*
x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x +
2*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*f*x
+ 2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 +
2*cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2*(co
s(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2*cos
(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*cos(2*
f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1))^2
)*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)), x) + 8*(a^2*c^2*f*cos(2*f*x + 2*e
)^2 + a^2*c^2*f*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*f*cos(2*f*x + 2*e) + a^2*c^2*f)*integrate((((cos(6*f*x + 6*e)*c
os(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e)
+ 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*cos(3/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e
))) + (cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x +
2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*sin(3/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*cos(5/2*arct
an2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - ((cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*
f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*cos(3/2*arctan2(sin(2*f*
x + 2*e), cos(2*f*x + 2*e))) - (cos(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(
2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)
*sin(3/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)
))/(((cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e)
+ 1)*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e
)^2 + 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*x + 6*
e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(2*f*x
+ 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2
*e)^2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*f*x +
2*e)^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 + 2*
cos(2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2*(cos(
2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2*cos(2
*f*x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*
x + 2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1))^2 +
(cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)
*cos(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e)^2
+ 2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*x + 6*e)^2
+ 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(2*f*x + 2*
e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^
2 + 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*f*x + 2*e)
^2 + cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 + 2*cos(
2*f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2*(cos(2*f*
x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x
+ 2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x +
2*e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1))^2)*(cos
(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)), x) - 2*(a^2*c^2*f*cos(2*f*x + 2*e)^2 +
a^2*c^2*f*sin(2*f*x + 2*e)^2 + 2*a^2*c^2*f*cos(2*f*x + 2*e) + a^2*c^2*f)*integrate((((cos(6*f*x + 6*e)*cos(2*f
*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x + 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*s
in(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*cos(1/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))) +
(cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x + 4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) -
2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*sin(1/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*cos(5/2*arctan2(si
n(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - ((cos(2*f*x + 2*e)*sin(6*f*x + 6*e) + 2*cos(2*f*x + 2*e)*sin(4*f*x +
4*e) - cos(6*f*x + 6*e)*sin(2*f*x + 2*e) - 2*cos(4*f*x + 4*e)*sin(2*f*x + 2*e))*cos(1/2*arctan2(sin(2*f*x + 2*
e), cos(2*f*x + 2*e))) - (cos(6*f*x + 6*e)*cos(2*f*x + 2*e) + 2*cos(4*f*x + 4*e)*cos(2*f*x + 2*e) + cos(2*f*x
+ 2*e)^2 + sin(6*f*x + 6*e)*sin(2*f*x + 2*e) + 2*sin(4*f*x + 4*e)*sin(2*f*x + 2*e) + sin(2*f*x + 2*e)^2)*sin(1
/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e))))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)))/(((
cos(2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*c
os(6*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e)^2 +
2*cos(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*x + 6*e)^2 +
4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(2*f*x + 2*e)
^2 + 2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2
+ 2*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*f*x + 2*e)^2
+ cos(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 + 2*cos(2*
f*x + 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2*(cos(2*f*x
+ 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x +
2*e) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*
e) + 1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*cos(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1))^2 + (cos(
2*f*x + 2*e)^4 + sin(2*f*x + 2*e)^4 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(6
*f*x + 6*e)^2 + 4*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e)^2 + 2*co
s(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(6*f*x + 6*e)^2 + 4*(
cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e)^2 + (2*cos(2*f*x + 2*e)^2 +
2*cos(2*f*x + 2*e) + 1)*sin(2*f*x + 2*e)^2 + 2*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 + 2*
(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*cos(4*f*x + 4*e) + 2*cos(2*f*x + 2*e)^2 + c
os(2*f*x + 2*e))*cos(6*f*x + 6*e) + 4*(cos(2*f*x + 2*e)^3 + cos(2*f*x + 2*e)*sin(2*f*x + 2*e)^2 + 2*cos(2*f*x
+ 2*e)^2 + cos(2*f*x + 2*e))*cos(4*f*x + 4*e) + cos(2*f*x + 2*e)^2 + 2*(sin(2*f*x + 2*e)^3 + 2*(cos(2*f*x + 2*
e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)*sin(4*f*x + 4*e) + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e
) + 1)*sin(2*f*x + 2*e))*sin(6*f*x + 6*e) + 4*(sin(2*f*x + 2*e)^3 + (cos(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) +
1)*sin(2*f*x + 2*e))*sin(4*f*x + 4*e))*sin(5/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1))^2)*(cos(2*f*x
+ 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e) + 1)^(1/4)), x))*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2
+ 2*cos(2*f*x + 2*e) + 1)^(3/4)*sqrt(a) - 16*(7*(5*a^2*c^2*sin(4*f*x + 4*e) + 4*a^2*c^2*sin(2*f*x + 2*e))*cos(
7/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)) - (35*a^2*c^2*cos(4*f*x + 4*e) + 28*a^2*c^2*cos(2*f*x + 2
*e) + 23*a^2*c^2)*sin(7/2*arctan2(sin(2*f*x + 2*e), cos(2*f*x + 2*e) + 1)))*sqrt(a))/((f*cos(2*f*x + 2*e)^2 +
f*sin(2*f*x + 2*e)^2 + 2*f*cos(2*f*x + 2*e) + f)*(cos(2*f*x + 2*e)^2 + sin(2*f*x + 2*e)^2 + 2*cos(2*f*x + 2*e)
+ 1)^(3/4))

Giac [F]

\[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\int { {\left (a \sec \left (f x + e\right ) + a\right )}^{\frac {5}{2}} {\left (c \sec \left (f x + e\right ) - c\right )}^{2} \,d x } \]

[In]

integrate((a+a*sec(f*x+e))^(5/2)*(c-c*sec(f*x+e))^2,x, algorithm="giac")

[Out]

sage0*x

Mupad [F(-1)]

Timed out. \[ \int (a+a \sec (e+f x))^{5/2} (c-c \sec (e+f x))^2 \, dx=\int {\left (a+\frac {a}{\cos \left (e+f\,x\right )}\right )}^{5/2}\,{\left (c-\frac {c}{\cos \left (e+f\,x\right )}\right )}^2 \,d x \]

[In]

int((a + a/cos(e + f*x))^(5/2)*(c - c/cos(e + f*x))^2,x)

[Out]

int((a + a/cos(e + f*x))^(5/2)*(c - c/cos(e + f*x))^2, x)